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Using the bound state version of the topological soliton model for the baryons we show that 
the existence of a bound (or quasi-bound) D-soliton state leads to the possibility of having hidden 
charm pentaquarks with quantum numbers and masses, which are compatible with those of the 
candidates recently reported by the LHCb experiment. The implications of heavy quark symmetry 
are elaborated. 
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I. INTRODUCTION 

The LHCb collaboration at CERN recently reported the discovery of two states, which in the quark model correspond 
to “pentaquarks” with hidden charm Ref. [if. The first one has the mass (width) 4380 ± 8 ± 29 (205 ± 18 ± 86 ) MeV 
and the second the mass (width) 4449.8 ± 1.7 ± 2.5 ( 39 ± 5 ± 19 ) MeV. The preferred quantum number assignments 
are J 71 ' = 3/2 _ ,5/2 + , respectively, although acceptable solutions are also found for additional cases with opposite 
parity. Note that their decay channel is J/if> p + which implies I = 1/2. After this announcement several articles that 
consider possible theoretical interpretations of the observed states have appeared. Some of these @41 suggest, with 
some variations, that the observed states may be interpreted as anticharmed meson (D or D*)-hyperon (E or E*) 
molecular states, while others base their description on diquark models |6H8]. The possibility that at least one of the 
observed peaks might be only a kinematical effect was discussed in RefA [9l4lll] . Further suggestions have been put 
forward in Refs. [13, 1131 . Moreover, the production and formation of the hidden-charm pentaquarks in 7 - nucleon 
collisions have been discussed in Refs. IldL 1151 . 

The purpose of the present note is to show that it is also possible to account for the quantum numbers and masses of 
these observed pentaquark candidates within the topological soliton picture of baryons. The possible existence of stable 
pentaquarks with negative charm in that framework was already considered many years ago fl6j . In this app roach the 
heavy flavor hyperons were described as bound states of heavy flavour mesons and a topological soliton [ 16H18I in the 
extension to heavy flavor of the bound state approximation (19h21) to the Skyrme model, supplemented with suitable 
symmetry breaking terms. This “Naive Skyrme Model” (NSM) formulation does however only approximately [22j 
incorporate heavy quark symmetry (HQS) according to which in the heavy quark limit the heavy pseudoscalar and 
vector fields become degenerate and, therefore should be treated on an equal footing [ 23 } . An improved way to proceed 
is therefore to apply the bound state approach to the heavy meson effective lagrangian [2ll-f27l , which simultaneously 
incorporates chiral symmetry and heavy quark symmetry (for details see Ref. [ 2 ^ and refs, therein). The possible 
existence of a C = —1 meson bound state in the context of a model consistent with HQS was first discussed in 
Refs. 0,(13]. There, however, only pseudoscalar mesons as described by the Skyrme model were considered in the 
light sector. In what follows we will refer to this as SMHQS formulation. It was later pointed out [32:] that the 
inclusion of light vector mesons in the corresponding effective lagrangian tends to push this state into the continuum. 
On other hand, the calculation in Ref. [33] . which incorporates the center of mass corrections in a more consistent way, 
still leaves the possibility of a loosely bound state. Additional arguments for the existence of C = — 1 meson-soliton 
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bound state have been given in Ref. [34j . 

In the NSM approach it is the Wess-Zumino term in the lagrangian, which is responsible for the difference in the 
interaction of the soliton and the mesons with opposite massive flavor quantum number. This term is repulsive in the 
case of mesons with massive antiflavor quantum number. In the case of S = +1 kaons, this rep ulsion, in combination 
with the repulsive effect of the meson kinetic energy term, pushes them into the continuum )35ll36|. (It should here be 
mentioned that the indications for the existence of the conjectured strange pentaquark “0(1540)” (37] have hitherto 
not been experimentally confirmed[38j). As the repulsive effect of the kinetic energy term weakens with increasing 
meson mass and the strength of the Wess-Zumino term is smaller for heavy flavors, the existence of anticharm (and a 
fortiori antibottom) meson-soliton bound states becomes possible. Below we show that the existence of a bound (or 
quasi-bound state) Z)-soliton state naturally leads to the possibility of having some hidden charm pentaquarks with 
quantum numbers and masses which are compatible with those of the candidates proposed in Ref.|l( 

II. BOUND STATE DESCRIPTION OF HIDDEN HEAVY FLAVOURED PENTAQUARKS 

Since the pentaquark candidates reported in Ref.QJ have no net charm quantum number we propose to describe 
them as bound states of a soliton and two pseudoscalar mesons (one charm and the other anticharm) in the present 
picture. We recall that in the bound state approximation while a bound C = +1 meson behaves as a quark [17H21| . 
a D -meson corresponds to a antiheavy-light quark pair [Ty, |3Q|, ElD] • To determine the possible quantum numbers and 
estimate the values of the associated masses we need to know which bound states of charm and anticharm mesons 
are possible. In fact, both the NSM and the SMHQS models lead to a number of meson bound states. Those bound 
states can be labelled by A; 71 ’, where k = i + £ and n = (— l) e+1 . Here, £ and i = 1/2 are the angular momentum and 
isospin of the bound pseudoscalar meson. In Table I the corresponding quantum numbers, binding energies b and 
hyperfine splittings constants c are given . The latter are needed for calculation of the nonadiabatic corrections to be 
discussed below. These results were obtained using the standard Skyrme model parameters /„• = 64.5 MeV, e = 5.45, 
which lead to the empirical masses for the nucleon and the A resonance. The values of the binding energies b and 
hyperfine splittings c associated with the C = ±1 mesons in the NSM scheme have been calculated in Refs. iTiii. fl8|. 
Those of the SMHQS scheme have been extracted from Ref.(29j| (Set 5) except for the binding of the C = — 1 state, 
which is taken from Ref.[33|. Note that most of the associated hyperfine splittings have not been given in those works, 
however. The results corresponding to the NSM have been obtained using the decay constant ratio fo/fir = 1.8. This 
value has been updated in recent years. The current estimate is fo/fir = 1-57 (JO]. While the use of this value hardly 
affects the predictions for the hyperfine splitting constants, one obtains an enhancement of about 100 MeV for all the 
C = +1 binding energies. Since the results for the binding energies obtained with fo/fir = 1-8 are closer to those of 
the SMHQS formulation we will assume that this overbinding is a consequence of the simplicity inherent to the NSM 
formulation. 


1 k n 

NSM 


SMHQS 

b(MeV) 

c 

b(MeV) c 

C = +l 1 1/2+ 

568 

0.20 

518 0.15 

0 1/2" 

355 

0.52 

239 0.30 

2 3/2 _ 

243 

0.15 

212 

1 3/2+ 

140 

0.28 

49 

1 1/2+ 

118 

0.03 

65 

C = -1 1 1/2+ 

38 

0.16 

54 


Table I: Quantum numbers of the meson bound states and associated binding energies and hyperfine splitting constants. 


The single meson spectra that follow from Table I are illustrated in Fig.l. It can be seen that the meson spectrum 
obtained in the NSM formulation is qualitatively similar to that obtained using the SMHQS formulation. Actually, 
this fact was already noticed long ago [22j |. 
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Figure 1: Binding energies and quantum numbers corresponding to possible soliton-(anti)charmed mesons. The left panel 
corresponds to the spectrum obtained using the NSM formulation while the right panel to the SMHQS one. The lowest 
= 3/2” configuration is indicated by the two black circles and the J w = 5/2 + by the open ones. 


Clearly, the C = — 1 meson has to be bound in the k"_~ = l/2 + state. Here, the subindex stands for the charm 
quantum number of the meson. For the C = +1 one, we have, however, several possibilities. The preferred quantum 
number assignments of the pentaquark candidates are J™ = 3/2 _ ,5/2 + , respectively. Since, obviously, 7 r = 7 r_ 7 r + it 
follows that 7 r + = —(+) for the 3/2 (5/2) pentaquarks. To determine k+ we note that, according to the usual rules, 
the total grand spin K = k+ + k- satisfies K = J + I where J and I are the total spin and isopin of the bound system. 
Noting that the observed states are isospin doublets, it follows that 1/2 < k+ < 5/2 for the state with J = 3/2 and 
3/2 < fc+ < 7/2 for that with J = 5/2. From Fig.l it is clear that the lowest lying meson configurations (fcT + ,fe^ _ ) 
that satisfy all the requirements are (l/2~, l/2 + ) for the J 77 = 3/2“ pentaquark candidate and (3/2 + , l/2 + ) for the 
J * = 5/2 + one. In the figure the first configuration is indicated by the two black circles and the second by the open 
ones. In principle many other states can be obtained by populating the alternative C = +1 bound states. The full 
list is given in Table II together with the corresponding masses up to non-adiabatic corrections. These have been 
obtained using 


M = M so i + 2m D (1) 

where M so i = 866 MeV is the soliton mass, mjj the pseudoscalar charm meson mass that we take to be rrip = 1867 
MeV and b± are the meson binding energies given in Table I. 

It is seen that three low lying 3/2“ with masses (4.21-4.32, 4.31-4.33) GeV for (NSM,SMHQS), respectively, are 
predicted in both schemes. The rather large width found in the experiment might be explained by these 3 close lying 
states. One the other hand only one single 5/2 + state at (4.42, 4.50) GeV is predicted. 

So far we have not included the non-adiabatic contributions to the masses. To first order perturbation theory the 
rotational corrections can be obtained by considering 

2 ^ < ( I(k + k-) K ) J \(R - 0) 2 | (I(k + k.) K ) J >, (2) 

where Q is the moment of inertia of the soliton and O the total isospin of the meson bound system. Moreover, R has 
the role of spin of the light quark system (basically the spin of the rotating soliton, which coincides with its isospin). 
For a system of a C = +1 and a C = — 1 bound mesons one has 0 = 0 + + 9-. The calculation of the matrix element 
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I 


k"~ 

K 

r 

NSM SMHQS 

1/2 

1/2+ 

1/2+ 

0 

1/2+ 

4.00 

4.03 




1 

1/2+ 






1 

3/2+ 



1/2 

1/2- 

1/2+ 

0 

1/2- 

4.21 

4.31 




1 

1/2- 






1 

3/2- 



1/2 

3/2" 

1/2+ 

1 

1/2- 

4.32 

4.33 




1 

3/2- 






2 

3/2- 






2 

5/2" 



1/2 

3/2+ 

1/2+ 

1 

1/2+ 

4.42 

4.50 




1 

3/2+ 






2 

3/2+ 






2 

5/2+ 



1/2 

1/2+ 

1/2+ 

0 

1/2+ 

4.44 

4.48 




1 

3/2+ 






1 

3/2+ 




Table II: Quantum numbers and 0(N°) masses (in GeV). 


in Eq.(© requires some assumptions and/or approximations. The usual procedure in the bound state approach (BSA) 
to the Skyrme model is to employ the approximation < Of >= c?£y(fc,; + 1). We denote this as BSA option. 

Here, Cj is the corresponding hyperfme splitting constant. As noted in Ref. [39[ this approximation does not hold in 
the HQS limit, where it can be shown that < Of >= 3/4, since for all the cases considered the bound mesons have 
isospin 1/2. We denote this as HQS option. Since it is not yet settled which the best way to proceed at the charm 
mass scale is, both cases will be considered here in order to determine the related uncertainty. The second one has to 
do with the way in which < R ■ 0 > is calculated (see Ref. 0). In fact, this term might induce mixings between states 
with different K and/or k+. This effect will be neglected in what follows. The resulting formula for the rotational 
corrections to the mass of a system composed by a soliton and two bound mesons (one in a state with and the 
other with fc_) is then 


Af ,Ji &+, i AT) — |d(A + 1) + c+c_ [K(K + 1) — fc+(&+ + 1) — k-(k- + 1)] + (5 

'c_|_ + c_ c_|_ — c_ fc+(fc+ + 1) — k-{k- + 1) 


[J{J + 1) - K(K + 1) - /(/ + 1)] 


K(K + 1) 


} (3) 


where 


cf_ fc+(fc+ + 1) + c 2 _ k-(k_ + 1) BSA 

S={ (4) 

1 3/2 HQS 

Using the BSA option one recovers the mass formula given in Refs. 0,0- Using the parameters of Table I the 
rotational corrections can then be calculated. Note that in the case of SMHQS scheme most of the hyperfme splittings 
needed for the calculation of the rotational corrections have not been given in the literature. Thus, only the predictions 
as obtained in the NSM scheme are reported in Table III. The quoted values correspond to the average between the 
results obtained using each of the two options for the calculation of the rotational corrections. The corresponding 
uncertainty is considered to be half of the difference between these two values, which turns out to be about 70 MeV. 
From this table we see that the prediction for the mass of the lowest 3/2 - is in the range 4.33 — 4.47 GeV. In the case 
of the 5/2+ the predicted mass is in the range 4.57 — 4.71 GeV. 
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r 

Mass [GeV] 

1/2+ 

4.11, 4.14, 4.51, 4.57, 4.59 

3/2+ 

4.16, 4.52, 4.60, 4.60 

1/2 

4.30, 4.35, 4.44 

3/2” 

4.40, 4.43, 4.48 

5/2~ 

4.50 

5/2+ 

4.64 


Table III: Masses including rotational corrections for the NSM scheme. All quoted values have an uncertainty of about 70 MeV 
due to the ambiguities in the formula for the rotational corrections. The lower value of this uncertainty range corresponds to 
the BSA option for the calculation of the rotational corrections while the upper to the HQS one. 


III. CONCLUSIONS 

The utility of the bound state interpretation described above is of course hinged on the existence of a C = — 1 
meson bound state. The calculations made within the framework of the NSM approximation indicate that such a 
state does exist. The SMHQS approach is less definite on this issue, but most of the calculations point towards the 
existence of a loosely bound state in the C = — 1 channel. It should be noted that even if it is unbound but lies close 
to threshold, one can still argue that an attractive DD interaction can make the whole soliton-Z1-.D system bound. 

Given the assumption of the existence of a loosely C = — 1 bound state, the existence of a pentaquark-type state 
with quantum numbers (I,J V ) = (l/2,3/2“) and mass in a region compatible with the LHCb observation follows 
naturally. Note that in the present picture at least one of the components of this state comes from populating the 
k+ + = 1/2“ meson state. In this sense it appears as akin to the A(1405). Note also that according to HQS such 
a meson bound state should be degenerate with a k^ + = 3/2“ state (see Fig.l). This probably then explains the 
existence of the other two 3/2“ state close by. The situation concerning the (J, J 7r ) = (l/2,5/2 + ) state is less 
clear. The model does predict such state but the associated mass lies about 4.6 GeV, which is somewhat too high as 
compared with the observed value. It has in fact been suggested that the observed peak at 4450 MeV might even be 
a kinematical effect [9j. In any case it is important to recall that, given the approximations made in the calculation 
of the masses, the quoted values have to be viewed only as first estimates. 

The present model obviously predicts the existence of several other states. In particular, it predicts two l/2 + and 
one 3/2 + states with masses in the range ~ 4.1 — 4.2 GeV, which arise by putting the C = +1 meson in the lowest 
k^ + = l/2 + bound state. Whether those states are too wide to be discriminated by the experiment remains to be 
seen. 

This bound state approach picture can be fairly straightforwardly extended to the bottom sector, where it would 
imply the existence also of hidden bottom pentaquarks in view of the large mass of the bottom mesons. This would 
probably require a full calculation within the SMHQS scheme, which is anyhow required to obtain more accurate 
predictions for the properties of the states discussed in the present work. 
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